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Abstract: Cross sections for top quarks provide very interesting physics opportunities, 

being both sensitive to new physics and also perturbatively tractable due to the large top 
quark mass. Rigorous factorization theorems for top cross sections can be derived in several 
kinematic scenarios, including the boosted regime in the peak region that we consider here. In 
the context of the corresponding factorization theorem for e'’'e“ collisions we extract the last 
missing ingredient that is needed to evaluate the cross section differential in the jet-mass at 
two-loop order, namely the matching coefficient at the scale fi czl mt. Our extraction also yields 
the final ingredients needed to carry out logarithmic resummation at next-to-next-to-leading 
logarithmic order (or N^LL if we ignore the missing 4-loop cusp anomalous dimension). This 
coefficient exhibits an amplitude level rapidity logarithm starting at 0(ag) due to virtual top 
quark loops, which we treat using rapidity renormalization group (RG) evolution. Interestingly, 
this rapidity RG evolution appears in the matching coefficient between two effective theories 
around the heavy quark mass scale /r ~ m*. 
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1 Introduction 


The top quark mass is one of the most important parameters in the Standard Model. As 
the heaviest observed fermion, the top quark provides an important probe for the Higgs sec¬ 
tor, and gives dominant contributions to many electroweak observables, thus providing strong 
benchmark constraints for extensions of the Standard Model. Furthermore, the mass of the 
top quark and the Higgs boson represent crucial parameters in studies of the stability of the 
Standard Model vacuum [1-6]. Precision measurements of the top quark mass are a difficult 
task due to challenges from both experimental and theoretical sides, mainly related to the fact 
that the top quark is a colored particle. 

The current value of the top quark mass from a combined analysis of Tevatron and LHC 
data is nit = 173.34 ± 0.76 GeV [7], see also [8, 9]. The precision obtained in this result relies 
on Monte Carlo (MC) based template and matrix element methods, which aim to account 
for essentially all of the kinematic final state information in the top quark events. However, 
this approach does not account for the relation of the extracted MC top quark parameter to 
an unambiguous field theoretic QCD top mass definition [10-12]. At the time of writing, no 
procedure to systematically quantify and improve this relation exists. While it seems unlikely 
that the template and matrix element analyses can be based on first principle QCD calculations 
which can be systematically improved to specify the top mass scheme unambiguously, it is quite 
plausible that other highly sensitive top mass observables can be devised which can clarify the 
issue by making high precision theoretical calculations feasible. 

One method to determine mt in a well-defined mass scheme from a kinematic spectrum 
with small uncertainties has been discussed in Refs. [10, 13, 14]. Here the hemisphere dijet 
invariant mass distribution in the peak region for the production of boosted tops in electron- 
positron annihilation was suggested as an observable and it was shown that hadron level pre¬ 
dictions of the double differential distribution can be carried out in a stable manner within a 
constrained set of top quark mass schemes. It was in particular demonstrated that the location 
of the peak of the distribution is highly sensitive to the top quark mass, and that only specific 
low-scale short-distance mass definitions are suitable for high-precision extractions. Although 
the effective theory setup developed therein was devised for the context of a future e"'“e“ col¬ 
lider, the approach can be extended to the environment at hadron colliders taking into account 
the complications related to initial state radiation, underlying event, parton distribution func¬ 
tions and dependence on jet algorithms and jet radius [15]. In Refs. [13, 14] the calculation for 
e'’'e“ annihilation was carried out at Next-to-Leading Logarithmic (NLL) accuracy with the 
perturbative ingredients at 0{as). In this paper we provide a result for the 0{a'^) matching 
correction at the scale ^ ^ mt for the e^e“-collider setup. Taken together with the known 
0{a‘^) results for the jet function in the heavy-quark limit from Ref. [16], for the massless 
soft function from Refs. [17-19], and input from previous form factor calculations for massless 
quark production [20, 21], our result provides the last missing ingredient needed to extend the 
e'’'e“ boosted top jet analysis to ©(o;^). In turn, with known results, these fixed order con¬ 
tributions can be accompanied with resummation of logarithms up to next-to-next-to-leading 
logarithmic order (NNLL). Up to the missing four loop cusp anomalous dimension, which is 
known to give a very small correction (see e.g. [22, 23]), all ingredients are also available for 
N^LL. 

Boosted top quark production with subsequent decays in the peak region of the invari- 
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Figure 1. Scales and effective theories with associated structures in the factorization theorem for 
boosted top production (Q ^ mt) with jet invariant masses close to the top mass. The superscripts 
(5) and (6) indicate the number of dynamic flavors in the theory. Note that in this context SCET just 
plays a role of an intermediate EFT with all invariant mass fluctuations above or of order the mass 
scale, in which the observable is not yet measured. For definiteness we also display the scaling of the 
EFT modes in light-cone coordinates. 


ant mass distributions involves physical effects in a range of widely separated energy scales. 
The hierarchy between the production energy Q, the top mass mt, the decay width T^ and the 
hadronization scale Aqcd is given hy Q ^ mt^Tt> Aqcd- Given this hierarchy of scales, the 
cross section contains large logarithms of ratios of these scales which spoil the perturbative ex¬ 
pansion in Us- This necessitates to replace fixed order computations by resummed calculations. 
The Effective Field Theory (EFT) setup devised in Ref. [13, 14] disentangles the fluctuations 
at the different scales and allows us to resum the logarithms through renormalization group 
evolution (RGE).^ 

We are interested in the peak region where each of the jet invariant masses, for the top st 
and antitop sj, is close to the top quark mass, i.e., 

Stt — rrii , ^ 

st,i^— - <.mt. ( 1 . 1 ) 

mt 

For this kinematic region both of the hierarchies Sti ^ Et and Sti ^ Ej are allowed. The 
sequence of the EFTs and the corresponding modes relevant for this problem are displayed 
in Fig. 1. First, hard modes with fluctuations with virtualities of order ~ Q are integrated 
out in QGD. The corresponding low-energy theory containing collinear and soft modes is Soft 
Gollinear Effective Theory (SGET) [26-29], which allows to resum large logarithms between 
Q and mt- In a second step all fluctuations with virtualities of order ~ mt are integrated out, 

^Note that our boosted top limit differs from the application of HQET in Ref. [24], which considers ti 
prodnction with slow top quarks. It also differs from the work of Ref. [25], which considers the top-pair 
invariant mass for boosted top quarks, rather than the individual boosted top jets. Hence the factorization 
theorem for our case differs from the ones considered there. 
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and SCET is thus matched onto boosted Heavy Quark Effective Theory (bHQET), an EFT 
with ultracollinear and ultrasoft modes at a lower invariant mass scale, which allows to resum 
logarithms between and The factorization theorem for the double differential cross 
section in e'’'e“ collisions reads 
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Here do denotes the tree level cross section for e+e" ^ qq. The terms Hq and Hm are hard 
functions related to the matching from QCD to SCET at the scale ^ ^ Q and from SCET 
to bHQET at the scale ix ^ m, respectively. The terms Jb and S denote the jet and soft 
functions, respectively, which are nonlocal matrix elements in bHQET. Note that we use Jb 
for the heavy-quark jet function, rather than the symbol B employed in Refs. [13, 14, 16]. 
Here Jb describes the dynamics of the ultracollinear radiation inside the t ox i jet at the scale 
H ^ St- The function S incorporates the ultrasoft cross talk between the two jets at the scale 
jjL ~ mst/Q, which is C>(Aqcd) in the peak region, and perturbative in the tail above the 
peak. In Eq. (1.2) the RGE between the characteristic scale of each function and the common 
renormalization scale are implicit. We stress that in SCET the top quark is considered as 
dynamical and hence the RGE takes place with six active flavors, while for the ingredients that 
arise in bHQET there are only five dynamical flavors in the evolution. Note that it is possible 
that the 0{mtaslQ) power corrections indicated in Eq. (1.2) are absent, but we are not aware 
of a rigorous proof at this time. 

It is through the residual mass term 5m appearing in the bHQET jet functions Jb that the 
top quark mass scheme is specified unambiguously beyond tree-level. For order-by-order stable 
perturbative behavior, the top quark mass scheme employed should be free of the O(Aqcd) 
renormalon ambiguity, thus excluding the pole mass (specified by 5m = 0) as a choice. Fur¬ 
thermore, the parametric scaling of higher order corrections defining the mass scheme must be 
set by scales associated to the measurement, namely ^mt, in order not to violate the 

power counting required for the factorization. This excludes employing the MS mass where 
these corrections scale as 5m ~ Oismt- Valid options include the jet mass scheme [13, 14, 16] or 
the MSR mass scheme [10, 16] which matches continuously onto MS. These two mass schemes 
have an adjustable cutoff parameter R which controls the scaling of higher order corrections. 

The exact algorithm to determine the two jet regions and the precise form of the observable 
is irrelevant for the structure of Eq. (1.2) as long as parametrically st ~ sj, but does matter 
for the explicit perturbative expressions of its ingredients. The restriction st ~ Sf avoids large 
logarithms of the form ln(st/sf), and is satisfied by variables designed to study the peak region 
of both jets, such as thrust. In the analysis of Ref. [14] all particles were assigned to either 
of the two top jets depending on which hemisphere with respect to the thrust axis they enter. 
Thus the observable considered was physically close to event-shape distributions. The analysis 
of Ref. [14] for this inclusive jet observable was carried out at NLL', i.e. including perturbative 
ingredients at 0(as) and NLL resummation. At the time of writing the hard function Hq, 
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the bHQET jet function and the soft function are already known up to 0{a‘^) [16, 17, 20] or 
beyond, while resummation can be carried out to N^LL.^ Thus, the only relevant correction 
missing to perform a N^LL analysis for the double hemisphere invariant mass distribution and 
similar observables in the peak region is the hard function Hm at 0 (q:^). This correction will 
affect the normalization of the differential cross section, while the shape of the cross section 
is determined mainly by the jet and soft functions. Here NNLL' refers to NNLL resummation 
with 0{a1) fixed-order matching and matrix element corrections. 

In this paper we carry out the computation of the 0{a'^) correction to Hm- In Sec. 2 we 
outline two methods to perform the computation. Instead of directly calculating the current 
matching factor between bHQET and SCET, we can also exploit the knowledge of the QCD 
heavy quark form factor calculated in Refs. [30, 31] and various properties of the EET to 
extract the hard function. In Sec. 3 we carry out the computation at 0{ag) using this method 
and show how to handle issues associated with the number of active quark flavors. This yields 
the result given in Eq. (3.8) in terms of the pole mass. In the two loop expression for Hm we 
find terms of the form 

“T.r.ln(g) (1.3) 

The large logarithm ln(Q^/m^) is induced by the separation in rapidity of soft mass-shell 
fluctuations with the scaling ~ {m,m,m) from collinear mass-shell fluctuations 

with ■,p~ ■,Pa_) ~ {m?IQ, Q, m). It can not be eliminated by a choice of p or summed by the 
RGE in p. This effect is directly related to virtual top quark loops which first appear at 0{a‘^), 
and has been discussed in detail in Refs. [32, 33] together with other subtleties concerning the 
incorporation of a massive quark in primary massless jet production in SCET. In Sec. 4 we 
will explicitly carry out the matching calculation for the 0{a^CFTp) correction with primary 
massive top quarks, and demonstrate how the amplitudes factorize into collinear and soft 
components which each involve a single rapidity scale. We show that this factorization is the 
same as that for massless external quarks, computed in Ref. [33], up to a different constant 
term that appears in the collinear corrections. The direct computation of the SCET soft and 
collinear diagrams at 0{a'^CFTF) can be performed elegantly by first computing the virtual 
correction for the radiation of a “massive gluon” at one-loop and performing in a second step a 
dispersion integral. In Sec. 5 we show how to resum the type of rapidity logarithm in Eq. (1.3) 
using the framework of the rapidity renormalization group established in Refs. [34, 35]. We 
also demonstrate that the residual scale dependence of Hm on p significantly decreases when 
employing the complete two-loop correction, and assess the impact of the rapidity logarithm. 
We conclude in Sec. 6. 

2 Setup and Notation 

As described in Refs. [13, 14] for the description of the peak region we first match QCD onto 
SCET, and then SCET onto bHQET. The relevant current operators needed to define the hard 

^So far the only missing ingredient for N^LL resummation (besides the four-loop cusp anomalous dimension) 
was the anomalous dimension of the jet function Jb, or equivalently of the bHQET current, which we have now 
extracted from a recent result in literature in appendix B. 
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functions in Eq. (1.2) are 


Jqcd ='ip{x)rf^l;{x), 

^SCET = XnSlTf SfiXn , 

JbHQET = V WnY^r>;YnWlK_ , (2.1) 

where r(( = 7 ^ and Eq = 7 ^ 75 . The jet helds Xn = Wn^n and Xn = describe the 

collinear radiation in SCET, and contain the massive collinear quarks and [36, 37] and 
Wilson lines Wn,n where in position space Wn{x) = Pexp(i 5 r J^dsn ■ An{ns + x)). The 
ultracollinear radiation in bHQET is described by the heavy quark fields hy_^ _ and by Wn,n- 
The wide-angle radiation in SCET is described by soft Wilson lines Sn,n, where in position 
space S'n(x) = Pexp (^ig f^dsn ■ As{ns + x)), and ultrasoft Wilson lines Yn^ are the analogs 
with ultrasoft gluon helds in bHQET. The difference between the SCET helds and bHQET 
helds is that SCET still contains soft and collinear huctuations at the top mass scale, i.e. 
the SCET helds contain mass mode huctuations which scale as ■,p~^Px) ~ and 

(Q, m^/Q, m) or (m?/Q,Q,m) which are absent in bHQET. This makes our EET above the 
top mass scale an SCETn type theory. There are six havors in the MS running coupling in 
QCD and SCET, and hve havors in bHQET. 

The notation above differs from Ref. [14] which used a hybrid of SCETj and SCETn, 
where the current operator was written as 

^SCET = XnY^SlTfSnYnXn • (2.2) 


Here the Wilson lines Sn,n describe exclusively soft mass mode huctuations and have ultrasoft 
zero-bin subtractions. In Eq. (2.1) the SCET operator only describes soft huctuations above 
and of order of the mass scale m, and not far below m. This simplihes the setup for the 
matching coefficient calculation, which in particular can be viewed as going from a six havor 
theory to a hve havor theory. 

The matching coefficients between these effective theories are dehned by 

xSceT [1 + 0(™/<3)] , (2.3) 

JsCe 7’ = cl”'* XbHQET [1 + 0(S/m)] . (2.4) 


Here both the currents and Wilson coefficients refer to the renormalized quantities. When we 
refer to the bare objects we will indicate this explicitly as e.g. in Eor all quantities 

we consider we use the renormalized coupling constant. When we want to separate the color 
structures of the matching coefficients we will do so in the following way; 




_j_ _j_ 





Oias) 


(2.5) 


In all the objects above the coupling is renormalized in the MS scheme with the number of 
dynamical havors, nj, being either ni or {ni + 1) as indicated by the superscript. Here ni is 
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the number of light quarks, and the additional flavor indicates the heavy quark (here the top 
quark). The choice for the number of flavors in each of the expressions above is motivated by 
the scales at which these objects live compared to the top mass. Note that we have kept the 
number of flavors appearing in Cm unspecified, as it can be expressed in either the ni- or the 
{ni + l)-flavor scheme. We will be explicit about which scheme we are using in the equations 
below. 

The hard functions in Eq. (1.2) are related to the Wilson coefficients via 

HQ{Q,fi) = \CQ\^ Hm(m,^,^?j =\Cm\^. ( 2 . 6 ) 


Here the dependence on Q in the hard function Hm appears due to the boost factor Qjm. 

In Eq. (1.2) all the functions live at their respective scales and are evolved to a common 
scale /Xfinai through renormalization group running. While the jet and the soft functions have 
convolution running [14], the large logarithms of the hard matching coefficients are summed 
by multiplicative evolution factors. 


HevoliQi fJ'&nal'i — Hq 


for Q, Hm — "m and //final < On the LHS the dependence on hq and Hm only comes 

from higher order corrections when the objects in Eq. (2.7) are truncated at a given order in 
resummed perturbation theory. The same is true for the rapidity scales vq and I'm, which are 
induced by the rapidity RGE that will be discussed further below and in Sec. 5.1. We will 
frequently drop these arguments that appear after the semicolon. The evolution factors here 
obey the RG equations 

an ym J \m / \m ) 

where is the anomalous dimension for the squared current in bHQET. 

Eqs. (2.3) and (2.4) suggest two different methods that one can use to calculate the 0{a‘^) 
piece of Cm or equivalently Hm- 


1) Indirect calculation using the known result for Cq and the matrix elements for the QCD 
and bHQET current operators in pure dimensional regularization: 

Using Eq. (2.3) and (2.4), and taking matrix elements of the operators with onshell top- 
quark states as in [13], we have 

/ fT{'n-l+^)\ _ ^{ni) j rjini) \ qn 

WQCD / ~ ^m WbHQET/ • 

Using the relation between bare and renormalized bHQET currents 

j \ _ r^ini) ,^(bare,ni)v /9 1 fli 

WbHQET/ ~ ^bHQET WbHQET / ’ 


-7 - 


we get 
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( 2 . 11 ) 


Note that the terms on the RHS involve objects with different flavor number schemes for 
the strong coupling, which must all be converted to n;-flavors to get ■ Here we work in 
dimensional regularization for both UV and IR divergences and renormalize the quantities in 
the MS scheme. With this regulator we can use the known two loop result for the heavy form 
factor (JTqcd) given in Refs. [30, 31]. The result for Cq is also known [20, 21] in MS, and the 
result for -^bHQET determined by RG consistency as discussed below. Loop graphs in 

bHQET factorize into ultrasoft and ultra-collinear contributions, and in general each involve 
at most a single dimensionful scale. The use of dimensional regularization for both the UV and 
IR, and employing onshell external quarks, imply that these loop corrections in bHQET are 
scaleless and vanish, such that (^bHQET = 1- general, the IR divergences in the QCD and 
bHQET matrix elements will match up, and the UV divergences in (^bHQET*^) eliminated 
by the counterterm -Z^^hqet- dimensional regularization with 1/eiR = l/euv; this implies a 
cancellation of 1/e poles between and -^bHQET- Thus we can use the simpler relation 



wqcd / 


^bHQET 


( 2 . 12 ) 


2) Direct calculation by matching the SCET and bHQET current operators: 

Using Eq. (2.4) we can also just directly compute the Wilson coefficient from a matching 
calculation, computing partonic matrix elements using the same IR regulator in SCET and 
bHQET, 





SCET / 


! y(ni) \ 
WbHQET/ 


-^SCET 

bHQET 


(2.13) 


These matrix elements are form factors in the respective theories which we denote by F. We 
will use the same notation for the color structures in the perturbative expansion of Escet and 
FbHQET as in Eq. (2.5). We define the relation between bare and renormalized SCET currents 
by 


! _ y{ni+l) / ^(bare,nj + l)> 

WSCET / ~ ^SCET WSCET / 


(2.14) 


As usual the bare currents are /r-independent, so from Eqs. (2.10), (2.13) and (2.14) the //-RG 
equation for can be written as 





7£'et T) - 7j,HQET 




(2.15) 


where the current anomalous dimensions are computed order-by-order from the counterterms 
in the standard fashion 


'SCET IV) Ml — M ^^^Il^SCET ’ 'bHQET — M j^^'^^bHQET ' 


(2.16) 








The anomalous dimension for the SCET current is known to 3-loop order [38]. Up to two loops 
the result reads 


i'SHtHQ. t^) = [- 4^-g + 6] + ( “%J*'‘* y{cM3-4.^ + 48C3] 




dvr 

+ CfCa 


r /268 

47r2^ 

\ . 961 

IItt^ 

l-(^- 


)Tq + ^ + 
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+ {ni + l)Ci?rp 


-m ^ 


(2.17) 


where Lq = ln[(—— iO)//r^]. The bHQET anomalous dimension can be derived using one 
of the consistency relations [14] for the factorization theorem in Eq. (1.2): 


Iv = 7bHQET + 7bHQET = + ^75 > 


(2.18) 


where 75 indicates the soft function anomalous dimension for one hemisphere. Using the results 
for 'yjg given in Eq. (41) of Ref. [16] and for 75 given in Eq. (19) of Ref. [39] (which can be 
derived via consistency from the two-loop jet function anomalous dimension [40]) we find 


7bHQET 
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dvr 


[— 4T -|- 4j -|- 




niCpTF 
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+ CfCa 




L + 


196 



+ C’(as)) 


(2.19) 


where L = ln[(—— i0)/m^]. Expanding the recently calculated anomalous dimension in 
HQET at 0{ag) [41, 42] we extract in appendix B also the three-loop coefficient, which has - 
to our knowledge - not yet been displayed in literature. 

As mentioned above, the two-loop expression of Cm contains large logarithms of the form 
ln(—m^/Q^) ~ 0{as) which cannot be resummed using the RGE in /r. They are 
rapidity logarithms and originate from a separation of the soft and collinear mass modes which 
have the same invariant mass but different rapidity. These rapidity logarithms only appear 
inside Hm, and not for the other soft, jet, and hard functions in Eq. (1.2). Our focus here 
will be on the leading rapidity logarithms, which start contributing with the 0{a‘^CFTF) 
piece. The latter comes from virtual top quark loops, and hence we only need to compute the 
correction while the bHQET graphs give no contribution for this color structure, 

bHQET ~ 

To set up the stage for rapidity resummation we can factorize the current operators and 
its matrix elements into products of soft and collinear diagrams. 


WSCET / ~ WSCET WSCET C \'-'sCET /"• ’ 

/ ^Ci+l) \ _ / ^{ni+1) \ j ^{ni+1) 1 / ^{ni+1) \ 

WbHQET/ ~ WbHQET/’^ WbHQET/s WbHQET/^ ’ 


( 2 . 20 ) 


where the {n, s, h} labels in bHQET indicate n-ucollinear, ultrasoft, and fi-ucollinear contri¬ 
butions respectively. Note that in order to split up these corrections we must choose an IR 
regulator which preserves the SCETn nature of the theory. We will regulate the IR divergences 
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using a gluon mass A, which thus differs from the use of pure dimensional regularization dis¬ 
cussed above for method 1. In SCETn the individual soft and collinear diagrams have rapidity 
divergences, and using the regulator of Refs. [34, 35] the coefficients will depend on a rapidity 
renormalization scale v. Thus Eq. (2.13) can be decomposed into individual contributions 
involving n-collinear, h-collinear, and soft amplitudes, 
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/ n-(ni + l) 
WSCET 




/ n-i'n-l) \ ’ 
WbHQET/* 


i = n,n, s . 


( 2 . 21 ) 


This leads to 




(jim) 


m, |U, — I cl?*} 




m 


'^m,n 



( 2 . 22 ) 


where we included the dependence on scales and renormalization parameters. Thus we see 
that the logarithmic dependence on the Q/m boost variable is factorized by the rapidity 
regularization parameter u into collinear factors that depend on Q and a soft factor which 
does not. To sum the rapidity logarithms we can follow the standard approach of matching 
and running. 


We define hard functions = 

functions obey related RG equations. 




The individual Wilson coefficient and hard 


d (n,) 


Cm 




h: 
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' \ IU, I J 


(2.23) 


The z/-anomalous dimensions appearing here can be computed directly from the SCET and 
bHQET counterterms and depend only on m and /x. Taking Eqs. (2.10) and (2.14) and intro¬ 
ducing individual counterterm factors for each of the collinear and soft component amplitudes, 
noting that the bare coefficients are ix-independent, and using Eq. (2.21) we get 


7, 


j7(m,^)=zzAlnCi"'] = 


'■^SCET / 
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Ain 1 1 

dv bHQET/ 


d, 

= zz—in 
av 


yirii + l) 
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zx—In 
dzz 


^bHQET,i ’ 


i = n,n, s . 


(2.24) 


As we will see in detail below, individual contributions on the right hand side of Eq. (2.24) 
contain IR divergences, but they will always cancel to leave an IR finite result for the 7 A, 
when we fully expand in either the rez-flavor or (n/ -|- l)-flavor scheme for the strong coupling. 


3 Two Loop Determination of Hm from QCD heavy form factor 


In this section we use the first method outlined in Sec. 2 to determine the bHQET matching 
coefficient, Cm at two loops. Erom Eq. (2.12) the ingredients we need are the UV renormalized 
QCD two-loop heavy quark form factor, (i7qcA^)’ ™ dimensional regularization and the SCET 
matching coefficient, In the following we abbreviate the appearing logarithms as 


L = In 




Lq = In 



(3.1) 
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From Refs. [30, 31] we extract the renormalized two loop QCD heavy quark form factor 
result in the high energy limit, evaluated at an arbitrary scale rn, abbreviating 
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(3.2) 


Note that we keep the 0{e) part of the one loop piece in Fqcd"*"^^ since it yields a contribution 
when considering the cross terms in the expansion of the ratio in Eq. (2.12). (One can avoid 

®Note that in Ref. [30] the counterterm for the renormalization of the coupling constant contains an extra 
factor r(l -I- e), so that also additional finite terms are subtracted compared to the conventional MS renormal¬ 
ization. 
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considering these cross terms and obtain the same answer by taking the logarithm of Eq. (2.12).) 
We remark that in these expressions the pole mass scheme has been used for the top quark 
mass m. 

The other ingredient we need is the well known two-loop expression for Cq, widely used 
in the SCET literature, and obtained with the aid of the massless form factor calculation of 
Refs. [20, 21], 
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(3.3) 


The remaining quantities in Eq. (2.12) are the coefficient we wish to determine, and 
the counterterm -^bHQET • contributions to these two quantities can be easily distinguished 

since -^bHQET terms with powers of 1 /e, whereas is given by the finite 0{e^) con¬ 

tribution. Therefore, it is straightforward to distinguish these two quantities unambiguously. 
Since we wish to determine these with n/ active flavors, we must convert the strong coupling 
in (77 qcd^^) to the n;-flavor scheme using the decoupling relation 

+ 0{al)\, (3.d) 


a 


^+i)(^) = a^)(/r) l + a^)(^) 


n(m^, 0) — 


a^)(/r)r^l 


37r 


where the one-loop vacuum polarization at zero momentum transfer for a massive quark pair 
is given by 


n(m^,0) = 


asiiLjTp ( 


37r 






r(e) = 


as{n)TF 

3f 


Lr 

e 


+ f + Y 2 ^ 


(3.5) 


We need to keep terms up to 0(e) in Eq. (3.d) since they contribute in the dimensional 
regularization scheme we are using when multiplying 0(as/e) IR divergent terms in Eq. (2.12). 
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Using these results in Eq. (2.12) we find the following expression for -^bHQET’ 
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This result can also be extracted from earlier literature using the consistency relation for RG 
running between and the soft and the jet functions in Eq. (1.2). In particular, the 1/e^ 
terms in Eq. (3.6) are given by a term involving the lowest order /^-function, and the square of 
the one-loop result (due to non-abelian exponentiation), while the 1 /e terms are directly related 
to the two-loop anomalous dimension given in Eq. (2.19). This provides a key cross-check for 
■^bHQET hence for our result below for Cm^\ 

After cancellation of the 1/e and 1/e^ terms in Eq. (2.12) with the help of -^bHQET’ 

remaining 0{e^) terms give the desired result for Cm^\ With the top-mass in the pole scheme 
we find 
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(3.7) 


Finally we arrive at the main result of this section - the result for Hm = |Gmr ia the n/-flavor 
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scheme with the top-mass in the pole scheme 
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(3.8) 


As anticipated, all of the logarithms in this expression are minimized for n ^ m, except for 
the contributions in the last line that involve the rapidity logarithm a^CpTp ln(Q^/m^). To 
understand the origin of this type of logarithm in the context of the renormalization group 
requires a further factorization of into soft and collinear pieces, as in Eq. (2.22). In the 

next section we will carry out an independent calculation of the 0{a‘^CFTp) terms in 77m 
This sets up the rapidity renormalization group analysis of this term, which can be found in 
Sec. 5.1. In Sec. 5.2 we present the result for with the top mass renormalized in the 

MS scheme. 


4 Direct Computation of the O^oICfTp) Result 
4.1 Ingredients for the Calcnlation 

In this section we perform a direct computation of the o^CfTf piece of the matching coefficient 
Cm{m,Q/m, fi) due to massive quark loops using the second method from Sec. 2. We carry 
out the calculation in analogy to Refs. [32, 33], where the corresponding contribution to the 
matching coefficient at the mass scale for massless external quarks (in the following called 
“primary”) was computed. In this section we extend the calculation to the case of primary 
massive quarks. 

Starting from Eq. (2.13) we note that for the a^CpTF massive quark term, the bHQET 
graphs expressed in the usual npflavor scheme do not give any contribution. The SCET graphs 
do contribute, and should be expressed in the same scheme for the strong coupling. Using the 
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decoupling relation in Eq. (3.4) we obtain in the notation of Eq. (2.5) 



j^iCpTp, ni+l) 
-^SCET 




(4.1) 


a. 


Ml 


In ({Q, m, A, /i) 


Stt 




.("i+i) , MO 


The second term on the right hand side accounts for the coupling conversion of the SCET 
form factor from (n; + 1) to ni flavors.^ As discussed in detail below, we will use a massive 
gluon as an IR regulator A, such that 0(e) terms in the coupling conversion in Eq. (3.4) can be 
dropped. For the remainder of this section we will drop the superscript (n; + 1) on the SCET 
form factors. 

We adopt the calculational method of Refs. [32, 33], where the two loop graphs containing 
a “secondary” massive quark bubble are calculated by starting with one-loop graphs describing 
the radiation of a massive gluon with mass M and applying in a second step dispersion relations 
to account for the gluon splitting into a pair of secondary massive quarks with masses m. The 
corresponding relation can be written as 


(_i)^w 

. Up„{m ,p 

+ ie 


(— 4 ) f — pm\ 

_ 1 r dM2 ^ ^ ~P^ J 

p"^ + ie TT J _ ]\^2 _|_ 


Im 


p^ + ie 


Tl(m^, 0). 


(4.2) 


Here Il{m^,p^) is the gluonic vacuum polarization due to the massive quark-antiquark bubble. 


("i^/) = -i{p‘^gi,„ - pf,p„)n{m‘^,p^)6^^ = / d'^x e*^’"^(0|rj;;^(x) J^(0)|0), (4.3) 


with the imaginary 


Im[n(m^,p^)] = 9{p^— Aw?) g'^Tp 


,2\-e 23-2'i7r(3-c^)/2 




d-2+ 


Aw? 

p2 


1 - 


4m ^ 

p2 


(4.4) 


We note that the same method can be applied to account for any kind of secondary particles 
by a corresponding choice of the polarization function H. Eq. (4.2) allows us to express the 
contribution to the SCET form factor due to the massive quark loops as 

- (n(m2,0) - ^i^^^iy(J,’^"“'^)(g,m,A), (4.5) 

^Note that the subscript —>• used here and elsewhere stands for the plain replacement of the 

couplings and does not involve any expansion based on Eq. (3.4). 
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where the “on-shell” form factor is 
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(4.6) 


In Eq. (4.5) A denotes the gluon mass acting as our IR regulator, which we distinguish from the 
gluon mass M used in the dispersion integration. Since total bare quantities are /r-independent, 
we do not add // as an argument to the components of bare quantities at a specific order. In 
the massive quark contributions to the coupling are renormalized with the onshell 
subtraction, i.e. is given in the scheme with n; dynamic flavors. In Eq. (4.5) the 

second term accounts for the change to re; -|- 1 dynamic flavors. The form factor itself is still 
unrenormalized, as indicated by the (bare) superscript. We perform the MS renormalization 
for the SCET current using Eq. (2.14). Incorporating Eqs. (4.5) and (2.14) into Eq. (4.1) the 
result for can be written as 


n{CpTp, ni) 
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p{l) 

-^SCET 




Here the 1-loop form factor is a UV and IR divergent amplitude, and -^gcET^^ 

SCET current counterterm in the (re; -|- l)-flavor scheme. Using the explicit form of n(m^,0) 
in Eq. (3.5) one can rewrite Eq. (4.7) as 



j^(OS,CpTp,ha.re) 

-^SCET 


(Q,m) -I- Z, 


(CpTp) 

SCET 


(Q,rre,/r) 


(4.8) 


-h ^n(rre^,0) - ^SCETe) > 

where we see explicitly that the dependence on the IR regulator is canceled. Note that we 
could have also carried out the computation employing the (re; -|- l)-flavor scheme to determine 
which involves converting the bHQET form factor from the ni to (re; -|- l)-flavor 
scheme. In this case the cancellation of IR divergences occurs in a different manner, and 
involves the 0{as) bHQET form factor. This approach is discussed in App. A. 

Note that nothing in Eq. (4.8) depends on the low energy bHQET theory. Therefore the 
result applies equally well to the case where one integrates out the heavy quark loop without 
approaching the jet invariant mass threshold st m? and matches onto a re;-flavor SCET 
theory instead of bHQET. In this case the matching coefficient only contains the contribution 
from the massive quark loop and receives corrections starting at 0{a‘^CFTp), so switching 
between the ni and (re; -|- I)-flavor schemes only affects the corrections at 0{ag) and beyond. 
This is in close analogy to the case of primary massless quarks discussed in detail in Refs. [32, 
33]. 

4.2 One-loop computation for secondary massive gluons 

Having laid out the basic framework in the previous section we now start with calculating the 
one loop SCET heavy quark form factors for a top-quark of mass m with a massive gluon of 
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Figure 2. Non-vanishing EFT diagrams for the computation of the hard current at 0{as ) with primary 
massive quarks and secondary massive gluons with masses m and M, respectively. Soft-bin subtractions 
are implied for the collinear diagrams. 


mass M to be used in the dispersion relation. The complete unrenormalized SCET result for 
the current form factor at 0(as) can be written as 

dcET”’(e.™.M) = . (4.9) 

'-V-' 

The correction with primary massless quarks has been already calculated in Refs. [32, 

35, 43-45] and reads in d = 4 — 2e dimensions 
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(4.10) 


where Lq = ln(— and Lm = ln('^). The corresponding one-loop counterterm in MS 




reads 


.a, 


^SCET 


dvr \ 


(4.11) 


Fig. 2 illustrates the SCET graphs with massive gluons needed to compute For the 

first three graphs in Fig. 2 the form factor contributions are defined as prefactors to the spinors, 
Fi^'^ Un,pl^Un,p' for i = n,n,s and are computed using the SCET Feynman rules for massive 
quarks given in Ref. [36]. 

Due to the eikonal structure the result for the soft diagram, is same as that for 

primary massless quarks [here = //^e'^^/(47r)], 


^a,bare) ^ _ 2ig‘^CFL^" 
For the n-collinear diagram we get 


d'^fc 1 1 1 
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{2fY [P _ Qf.+ _ r^j^- [A;- -k ie] [k^ - M'^ + ie] ' 


(4.12) 


(4.13) 


We can decompose this contribution into a correction corresponding to the diagram with 
primary massless quarks, and a UV and IR-finite difference of terms which can be computed 
in 4 dimensions. 


p{l,bare) _ p{l,bare) /'pi(l,bare) _ p(l,bare)\ 
^n,m=0 ' \-^n ^n,m=0 ) ' 


(4.14) 
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After performing a contour integration in /c"*", carrying out the /c^-integration and rescaling 
the label momentum as k~ = zQ, the finite correction due to the mass of the primary quark 
yields 
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In SCET loop graphs include soft 0-bin subtractions [46] which ensure that there is no double 
counting of infrared regions. For the soft 0-bin subtraction of the dependence on the 

primary quark mass drops out, and we obtain the same result as for primary massless quarks, 
which is therefore fully contained in Note that the result in Eq. (4.15) does not 

contain any rapidity divergences, so that rapidity logarithms arise only in the computation of 
This can be understood from the fact that the corrections due to soft modes are the 
same for massless and massive primary quarks, so that the rapidity divergences in the soft 
sector and, by consistency, also in the collinear sectors have to agree in both cases. 

The n-collinear diagram corresponds to switching k~ and k'^ in Eq. (4.13). We perform a 
decomposition analogous to Eq. (4.14), 

— ^n,m=0 + Wn ^n,m=0 )■ l^.l/j 


The difference correction due to the primary quark mass is again UV and IR-finite and does 
not contain any rapidity divergences. Thus it yields for any choice of regulator the same result 
as the n-collinear correction, i.e. 

p(l,bare) _ p(l,bare) _ p(l,bare) _ p(l,bare) / . , 

^n,m=0 ^ n,m=0 ' 


Finally, we also have to consider the wave function corrections. In analogy to the computation 
in Ref. [14] we have 

w(i)_o- 2 ^ ~ 2 ei f Qm'^{3-e)-{Q'^k++ Qp'^ + m‘^k-){l-e) 

ztg 2 J (2TrY Q^[k^ - M2 + ie] [{k + pf - rrF + ie] ' ^ ’ 

Usingand decomposing the integrals into elementary one- and two-point functions 
we obtain 
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(4.20) 
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which uses the loop integrals 
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(4.21) 


The wave function renormalization constant ' is defined by taking the on-shell limit A —?■ 0 
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'2 Q 
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(4.22) 


where is the one-loop renormalization constant for the quark mass m in the pole mass 

scheme for the interaction with a massive gluon (with mass M). The wavefunction correction 
Z^^'^ can be written in terms of the wavefunction correction for primary massless quarks and a 
UV and IR finite remainder, 
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(4.23) 


The remainder contribution in d = 
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4 dimensions reads 
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- 140^ + 3a^) + 0(e) , 


(4.24) 


where a was given above in Eq. (4.16). 

The complete finite correction at one-loop, which accounts for the mass of the primary 
quark is given by the sum of the terms from Eqs. (4.15) and (4.24), 

M) = 2 (m, M) - {zf - (m, M). (4.25) 

This result will be used for our two-loop computation in the next section. 


4.3 Two-loop computation for secondary massive quarks 

In this section we use the one-loop results from Sec. 4.2 to calculate the two-loop graph with 
the massive quark loop, and to determine the CpTp contribution to Cm- First we compute 

using the one-loop result in Eq. (4.9). Again we can decompose 
the two loop SCET form factor into a primary massless component and a correction for primary 
massive top quarks: 

pfOSjCpT^jbare) _ p(OS,Ci?T^,bare) X p{CpTp) /a r\f-\ 

-^SCET “■^SCET,m=0 
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The calculation for primary massless quarks has already been performed in Ref. [33]. We 
display the result here for convenience; 
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(4.27) 


The contribution from the two-loop MS counterterm is known from the massless quark case 
and reads 
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(4.28) 


where L and Lq are defined in Eq. (3.1). The l/e'^ divergences in Eqs. (4.27) and (4.28) differ, 
and are reconciled only once we account for the additional scheme change correction in the 
last term of Eq. (4.5). The 6Fm ^ ^ term can be computed by inserting the one-loop massive 
gluon correction term of Eq. (4.25) into the dispersive integral (4.6) which can be performed 
in four dimensions. The result reads 


™ V 47r y \ 81 27 3 / 


(4.29) 


Thus the only modification in the massive quark loop contributions to the form factor for 
primary massive quarks with respect to primary massless quarks is a simple constant term. In 
particular no additional rapidity logarithm ~ ln(Q^/m^) appears, which can be again traced 
back to the universality of the soft corrections for massless and massive primary quarks. 

Assembling all the pieces above in Eq. (4.8) we get the following result for ■ 



which matches exactly with the CpTp result we obtained above in Eq. (3.7). In the next 
section we decompose the SCET form factor result into soft and collinear pieces in order to 
find the terms needed for the rapidity RGE analysis. 


4.4 Two Loop Ingredients for the Rapidity Renormalization Group 

In order to determine the ingredients needed for the rapidity renormalization group analysis, 
we now calculate the 0{a‘^CpTp) SCET form factor contributions for the individual collinear 
and soft sectors using dispersion relations. We will employ the symmetric ry-regulator [34, 35] 
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to regulate the rapidity divergences in the individual sectors. This corresponds to modifying 
the Wilson lines in the respective sectors according to 
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(4.31) 


and similarly for Wn and Sn- Here denotes the label momentum operator and w{v) is a 
dimensionless book keeping coupling parameter satisfying 
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(4.32) 


The one-loop form factor corrections for the radiation of a massive gluon have been already 
calculated in Ref. [35] for massless quarks. Including the modihcation due to the quark mass 
in Eq. (4.25) they read after expanding in rj 
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In the collinear results we have included the wave function contributions and Z^-/2. 

The soft-bin subtractions in the collinear diagrams vanish for the 77 -regulator. 

In direct analogy to Eq. (4.5) the corresponding two-loop expressions for the individual 
soft and collinear sectors read 
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(4.34) 


for i = n,n, s. Note that for this relation to make sense also the one-loop form factor corrections 
with a massless gluon have to be decomposed according to Eq. (2.21). To achieve this goal 
we use a gluon mass A <C as an infrared regulator which allows us to use the results in 
Eq. (4.33). As discussed in Sec. 2 , we absorb all divergences of the form 1 / 77 , rf /e^ in the form 
factors into separate counterterms -^gcET^i each sector, so that 
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The explicit results for the counterterms at one-loop are given by^ 
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while at two-loop they read 
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Note that the snm -^scET^n + -^scET^i + “^scet^s reproduces the result for the SCET current 

(C^ T 'l 

counterterm ™ (‘1-28). These results for the individual collinear and soft conn- 

terterms provide the necessary ingredients for determining the rapidity RGE for the collinear 
and soft sectors below in Sec. 5.1. 


5 Rapidity Evolution and Numerical Results 
5.1 Rapidity Renormalization Group Evolution 

In our result for the matching coefficient between bHQET and SCET at O(a^), given above 
in Eq. (3.8), we encountered a large logarithm a^GirTp ln(m^/Q^). We discussed the setup 
for the resummation of such logarithms above in Sec. 2 . As shown in Sec. 4 these rapidity 
logarithms are only related to contributions of the virtual massive quarks that appear in 
the glnon vacunm polarization, and hence are the same as in the threshold corrections for 
massless primary quarks in Ref. [33]. There it was anticipated that they can be resummed by 
exponentiation, as is common for these kinds of logarithms. Eor example, for the radiation 
of a massive gauge boson the rapidity renormalization gronp implies that this exponentiation 

® Although the full e-dependence in the expression proportional to 1/j; should be in principle kept unexpanded, 
this is only relevant to ensure that the coefficient of the l /?7 pole is explicitly /r-independent, which is also true 
order by order in its e expansion. Therefore we show here only the terms up to 0{e^) which contain the 
information we need later for the anomalous dimensions. 
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occurs to all orders in perturbation theory [32, 35, 43, 44], The difference in our case is that 
the rapidity logarithms start at two-loops, and hence involve the additional issue of one-loop 
induced corrections due to the scheme change in the coupling constant. 

Here we will show explicitly how to treat the rapidity logarithms at 0{a‘^CFTp) in a 
rapidity renormalization group framework, and subsequently demonstrate that they indeed 
exponentiate. We start from Eq. (2.24). Up to 0{a1) we only have a contribution from the 
CpTp dependent terms, 

= z^^lnZscET.i - z^^ln ZbHQET,i 


, d (CpTf) 
'^dj/ SCET,i 


(ni) 




where the second term accounts for coupling conversion from the (n/ -|- l)-flavor to n/-flavor 
scheme. As before, in the n/-flavor scheme the bHQET graphs give no contribution. The results 
from Sec. 4.4 can now be used to compute this i/-anomalous dimension. Using Eq. (4.36) we 
can calculate the one-loop correction, 
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which exhibits dependence on the infrared gluon-mass regulator A. The two-loop term above 
can be calculated using Eq. (4.37) which gives 
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where is defined in Eq. (3.1). Together these results determine the zy-anomalous dimensions: 
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Note that the IR regulator has canceled out, and that here the coupling [q:s(/u)]^ can be taken 
in either the ni or {ni -|- l)-flavor scheme since the anomalous dimension starts at 0{a'^) and 
the difference is higher order. This result suffices for solving the Z/-RGE equations at NNLL 
order. Using Eq. (2.22) and Eq. (2.23) we can write an analog of Eq. (2.7) for the zy-evolution 
of Hm- From Eq. (2.22) we have 
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With rapidity evolution this becomes 
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where on the LHS the dependence on i'q and r'm comes from truncating the resummed pertur¬ 
bation theory for objects on the RHS. Here the functions Hm^n = are given 
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and contain no large logarithms for ^ ~ m, and for i/g ~ Q and Vm — m, respectively. The 
evolution factor Hrrg sums the rapidity logs between Vm and vq, and is defined as follows 
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The general result for Vrrg, and the result at NNLL, will be given below. 

Similarly to the anomalous dimensions, we can also determine individual ^u-anomalous 
dimensions for the collinear and soft sectors, i = n,s,n, 


I fl, I 


= h^ln^SCET,i - +^lnZbHQET,i 


(5.10) 


- 24 - 



Repeating the steps below Eq. (5.1) we find 
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whose sum yields the same result for the 0(0^0fTf) /i-anomalous dimension of as the 
difference of Eqs. (2.19) and (2.17), 
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with Lm and Lq defined in Eq. (3.1). 

Eqs. (2.7) and (5.6) together include the evolution connected to Hm in the 2-dimensional 
plane, including that from invariant mass scales fim to fig, that from invariant mass 
scales fim to fi^nah and that from rapidity scales fq to Fm- As demonstrated in Ref. [35] the 
combined fi-F evolution can be performed along any path and the path independence implies 
the consistency equation: 
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However, similar to the example of the massive Sudakov form factor considered in Ref. [35] 
we can see from Eq. (5.4) that contains potentially large logarithms ln(/i/m) for an 
arbitrary path in //-z^-space. This is resolved by a prior resummation exploiting the fact that 
the derivatives in Eq. (5.13) are proportional to the cusp anomalous dimension. Since Cm is a 
matching coefficient between a {ni + l)-flavor and n;-flavor theory, we can express Eq. (5.13) 
in terms of the difference between the cusp anomalous dimensions rcusp[as] in the (ni + 1 ) and 
n;-flavor schemes. So for 7 ^^^ we obtain 
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which can be checked using the explicit perturbative expression of rcusp[as] up to two 
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Integrating Eq. (5.14) in fj, we obtain the resummed result for 
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Here the integration constant is the correction in the anomalous dimension 7 ^ 2 ^ 

that does not multiply a logarithm ln(|U^/m^). We are now in the position to write down a 
general expression for VaRc- Using Eq. (5.9) we find the all orders result 
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At NNLL order with the counting as(/r)ln(r'm/z^Q) ~ 1; we can expand this exponential to the 
first non-trivial order. At the order we are working 
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as can be seen from Eq. (5.4), where we have for definiteness employed the (n; + l)-flavor 
scheme. The evolution function u at NNLL accuracy reads 
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where r = /a^J^^\no) and the coefficients /3j and Lj are evaluated with nj flavors. 

To extend the analysis to N^LL resummation, one needs the result for the z^-anomalous 
dimension at 0{a^), which can be inferred from the coefficient of the rapidity 

logarithm appearing in a related DIS calculation [47] due to consistency (see Ref. [48]). 


5.2 Numerical Results 

In this section we explore the impact of the two-loop correction to the hard function Hm on 
the differential cross section and the corresponding improvement to the perturbative uncer¬ 
tainties. To do this we examine the evolved hard function i7evoi(Q) "ij/^finai! MQjt'm) 
from Eq. (2.7). This function fully captures the multiplicative contributions for the differential 
cross section factorization theorem in Eq. (1.2), including the matching at /ig ~ Q in 

the RG evolution from /rg down to Hm — m in the matching at encoded in Hm, 

and through the RG evolution from down to a scale /Ugnai where the soft and jet 

functions are evaluated.® Since the ingredient that has not been previously analyzed is Hm we 
focus our numerical study on the impact of this function and the associated reduction in the 
resulting fim dependence. For Hm'^\rn,Q/m, we employ Eq. (5.6), which provides 

®The soft or jet functions also contain an additional evolution which is not purely multiplicative [13]. This 
evolution affects the shape of the du/dsidst distribution and was evaluated up to NNLL' order in Ref. [16]. 
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a decomposition of this function into collinear and soft components, with i = n,n, s, plus 

a kernel I^irg which carries out the RG evolution in rapidity from i/q ~ Q to — w,. 

We begin by converting the result for the collinear and soft components in Eqs. (5.7) 
and (5.8) from the pole-mass scheme to the MS mass scheme with ni + 1 dynamic flavors via 

mpole = m("'+')(//) (l - ^ ^ ^5 20) 

The MS scheme is an appropriate renormalon-free short distance mass scheme to be employed 
in the hard function Hm- For consistency we also convert the results in Eqs. (5.7) and (5.8) to 
the {ni + l)-flavor scheme for the strong coupling. Together this yields up to 0{ag) 
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where Em = ln(?fi^//i^) and m = m("''’'^)(/i) is the MS mass for re; + 1 active flavors. For the 
bHQET evolution function when using the MS mass scheme, we expand the pole mass 

appearing in the anomalous dimension in Eq. (2.19) in terms of fht{fht) to obtain 
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where L = ln[(—— iO)/fh‘^]. For the z^-anomalous dimensions the MS results are obtained 
by the simple replacement m ^ fh, since they start at two-loops. For our central results below 
we use Hm = Vm = rnt and fig = i'q = Q. 

For our numerical analysis of Ffevoi we employ scale choices that are appropriate to the 
peak region of the differential cross section within bHQET. We hx Q = fJ-Q = 1 TeV, which 
is a possible c.m. energy for a future linear collider, and /ifinai = 5 GeV corresponding to 
the scale of the soft radiation. We do not vary these two scales here since their impact 
and associated uncertainties have been analyzed elsewhere [14]. They matter only for the 
overall normalization and thus cancel in the normalized spectrum. In addition we use the MS 
mass rht{mt) = 163GeV or pole mass mt = 171.8 GeV using the two-loop conversion, and 
oi^s\'mz) = 0.114 [23, 49] and using two-loop conversion aX ^ = rht to obtain For 

results with RG evolution that sums large logarithms we use the so called primed counting, 
i.e. our results at NLL' and NNLL' include NLL and NNLL evolution kernels together with 
the hard function boundary conditions at 0(as) and C>(a^), respectively. ^ For the rapidity 
evolution we use the expression in Eq. (5.17), and the default rapidity scales vq = Q and 
Vm = mt, where mt is either the MS mass mt{mt) or the pole mass. 

To determine the impact on the normalization we first note that the two-loop fixed order 
corrections to turn out to be small, giving at the central scale = mt{mt) a 2% 

correction and the hxed-order series 

(mt, = mt] = 1 -|- 0.126(l-loop) -|- 0.015(2-loop) = 1.141. (5-24) 

\ mt / 

In the top-left panel of Fig. 3 we display the evolved hard function 77evoi at the hrst three 
orders in resummed perturbation theory for values of in the range fht/2 < < 2fht. We 

use the MS mass scheme and the expressions for 77^*^^\ and from Eqs. (5.21) 

and (5.22). As already observed in Ref. [14], there is a significant correction when going from 
LL to NLL' order which more than doubles 77evoi- From NLL' to NNLL' we observe that the 
correction is notably smaller, indicating that the series has stabilized. Although the magnitude 
of these corrections is not captured by the /im variation, it is of the size expected from studying 
the uncertainty associated to the /ignai variation. The complete study of the /ifinai variation 
requires including the jet and soft functions, which cancel the /ifinai dependence of 77evoi to the 
order one is working. We leave this for future work rather than taking it up here. We observe 
that the Hm dependence signihcantly decreases as we go to higher order. This behavior is shown 
best in the top-right panel of Fig. 3, where the same curves are plotted, but now normalized 
to 77evoi(7im = mt) at the respective order. The two-loop result for the hard function 
plays a key role in this reduction of the scale dependence at NNLL'. Note that the size of the 
^rn variation of the blue dashed curve at 2% correlates well with the size of the NNLO fixed 
order correction in Eq. (5.24), which gives a +2% correction. Therefore it is reasonable to take 
the /Um variation of the solid red curve in this hgure as an estimate of the 0{ag) correction in 
Eq. (5.24), which we take to be ±0.2%. 

^Going from NNLL' to an even higher order in the resummation, N^LL, does not affect any of the conclusions 
in this section, and therefore, for convenience, we carry out our numerical analysis at NNLL'. 
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Figure 3. Upper panels: Plots of the residual dependence on the matching scale /im for the unnormal¬ 
ized (left) and normalized (right) evolved hard function TJevoi at three different orders in the evolution, 
using the MS mass. Lower left panel: Comparison of the scale dependence at NNLL' for the MS mass 
and the pole mass. Lower right panel: Impact of varying the ratio of rapidity scales vqjvm by a factor 
of two at NNLL' as a function of Hm, with the MS mass. 


In the lower-left panel of Fig. 3 we compare the dependence on at NNLL' for the 
MS mass with the corresponding result for the pole mass. In the pole mass case we em¬ 
ploy Eqs. (5.7) and (5.8) for Hm,n, and . We see that the pole mass exhibits a 

larger sensitivity to the renormalization scale implying a slightly slower convergence of the 
perturbative series, potentially related to IR renormalon effects. 

Finally, we can analyze the impact of the terms related to rapidity logarithms. For = 
fht{rht), these terms yield a numerical contribution of —0.0014 in the fixed-order full hard 
function fim = hfi) in Eq. (5.24). Due to a relatively small coefficient, they 

do not give a significant correction in comparison with the remaining two-loop contributions 
which give a numerical correction of 0.0166. Therefore, we anticipate the dependence on the 
rapidity scales i^q and Vm to be rather mild. In the lower-right panel of Fig. 3 we plot iLevoi at 
NNLL' for the MS mass as a function of fj-m, but now with three choices for vqjyrn- To obtain 
these results we varied vq up and down by a factor of two, but we note that equivalent results 
are obtained by instead varying Vm by a factor of two. We see that varying uq jvmhy a. factor of 
2 gives a negligible effect compared to the residual dependence at this order. Therefore, we 
conclude that including an uncertainty from zz-variation is not necessary to obtain an estimate 
of the overall perturbative uncertainty of the cross section. 
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6 Conclusions 


In the context of EFT factorization for boosted top quark production, we have extracted the 
hard function Hm = jCmP describing virtual fluctuations at the top-mass scale, completely 
at two-loop order using earlier results from Refs. [30, 31]. This result provides the last miss¬ 
ing ingredient needed to make N^LL resummed predictions (up to the 4-loop cusp anomalous 
dimension) for the invariant mass distribution of top-jets in the peak region using the factoriza¬ 
tion theorem of Refs. [13, 14] given in Eq. (1.2). Particular focus was given to the contributions 
to Hm from heavy quark loops, which induce terms with a large logarithm a^CFTp\'a.{Q‘^/m?') 
that can not be treated with standard RG evolution in These terms were computed once 
more directly using collinear and soft matrix elements in SCET, and we have shown how they 
can be factorized using a rapidity cutoff v, and RG evolved using rapidity renormalization 
group equations. Interestingly, this factorization and RG evolution occurs within the Wilson 
coefficient Cm and hence at the amplitude level. Using our result for Hm we have assessed the 
remaining perturbative uncertainty associated to the top-mass scale, ^m — and estimate it 
to be very small, ±0.2%, predicting that the two-loop result for Hm provides a very accurate 
result for this function. The total normalization uncertainty in the differential cross section is 
expected to now be dominated by that from 0(af) perturbative corrections to the low-scale 
soft and jet functions, which could be estimated by a dedicated study of the residual /Ugnai 
dependence at NNLL' order. 
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A Direct Calculation of Cm in the (n/ ± l)-flavor scheme 


In Sec. 4 we directly computed the 0{a'^CFTF) massive quark correction to Cm’’'^ by using 
form factors in the n;-flavor scheme. Since this coefficient lives at the border between the 
{ni ± 1) and n;-flavor theories, we could just as well have carried out the calculation for Cm by 
using form factors in the (n/ ± l)-flavor scheme, and then converted to an n;-flavor coupling 
at the very end. Of course the same result is obtained in this approach, but there are a few 
subtle differences in the calculation, which we discuss here. 

In particular, in Sec. 4 we noted that for the 0{a‘^CFTF) correction in the n;-flavor scheme, 
the bHQET graphs give no contribution. However, using the (ni ± l)-flavor scheme for the 
strong coupling this is no longer the case. To see this, consider the ratio in Eq. (2.13) and 
express the denominator in the (n; ± l)-flavor scheme by inverting the decoupling relation given 
in Eq. (3.4): 
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Expanding in ag and using the notation in Eq. (2.5) we then get 
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Here the second term comes from converting the strong conpling constant to (n/ + l)-flavors in 
the one-loop bHQET graph. Below we drop the flavors snperscript on the form factors. Here 
it shonld be understood that all the terms are now expressed in the (n/ -|- l)-flavor scheme. 
Then combining Eq. (A.2) and Eq. (4.5), and Eq. (2.14) we get 
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Note that both and are IR divergent. This result can be simplified by noting 

that in any flavor scheme the one-loop Cm'^ is given by the difference of one-loop renormalized 
SCET and bHQET amplitudes: 

,/x^ = Eg(]grp(Q,m,A,^) — E^jjQp,rp^— ,A,. (A. 4) 

Using Eq. (A.4) in Eq. (A.3) we can then write down a simpler expression for 
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This result can be used to compute FjFFTF,ni+i)^ Comparing it with Eq. (4.8) we see that it 
can be rewritten as 


CjCrT,. „,+i) ^ ci» , (A,6) 
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and hence is fully consistent with determining from Eq. (4.8) and then simply 

applying the coupling conversion in Eq. (A.l) in the result. Note that in this (n; + l)-flavor 
scheme approach the bHQET one-loop amplitude contributes and plays an important role in 
obtaining the scheme conversion term involving in the last line of Eq. (A. 6 ). 


B bHQET current anomalous dimension at 0{a^) 


To extent the resummation of large logarithms in the factorization theorem in Eq. (1.2) from 
NNLL to N^LL the only missing ingredient - besides the cusp anomalous dimension at four- 
loops - is the 0{a^) noncusp anomalous dimension of the bHQET jet function or equivalently 
of the bHQET current (which are related to each other via Eq. (2.18) with the known three 
loop result for 75 ). The latter has not been so far given in the literature, but can be extracted 
from a recent result for the three-loop anomalous dimension of a cusped Wilson loop [41, 42], 
which is equivalent to the full anomalous dimension in HQET. Expanding their result in the 
lightlike limit x ~ m/Q —>■ 0, we obtain with the help of the Mathematica package HPL [50] 
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where L = ln[(—— z0)/m^]. The coefficient of this logarithm is proportional to the well- 
known lightlike cusp anomalous dimension at three loops, Tcu^p, while the non-logarithmic 
ingredient of Eq. (B.l) represents the noncusp part. Together with the corresponding anoma¬ 
lous dimension of the SCET current this enables one to predict the logarithmic structure of 
Hm at three loops by solving Eq. (2.15). Eurthermore it allows one to extract the last missing 
ingredient to predict the full IR-divergent structure of the three-loop full QCD form factor for 
massive quarks for m <C Q, which is for example in Ref. [31] the coefficient in Eq. (63). 
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